home *** CD-ROM | disk | FTP | other *** search
/ Sprite 1984 - 1993 / Sprite 1984 - 1993.iso / src / lib / c / etc / isnan.c < prev    next >
C/C++ Source or Header  |  1991-08-30  |  2KB  |  112 lines
  1. /* 
  2.  * isnan.c --
  3.  *
  4.  *    Machine-dependent procedure to determine whether a double is a
  5.  *    "NaN" floating-point number.
  6.  *
  7.  * Copyright 1989 Regents of the University of California
  8.  * Permission to use, copy, modify, and distribute this
  9.  * software and its documentation for any purpose and without
  10.  * fee is hereby granted, provided that the above copyright
  11.  * notice appear in all copies.  The University of California
  12.  * makes no representations about the suitability of this
  13.  * software for any purpose.  It is provided "as is" without
  14.  * express or implied warranty.
  15.  */
  16.  
  17. #ifndef lint
  18. static char rcsid[] = "$Header: /sprite/src/lib/c/etc/RCS/isnan.c,v 1.2 91/08/29 21:52:22 rab Exp $ SPRITE (Berkeley)";
  19. #endif /* not lint */
  20.  
  21. #include <math.h>
  22. #include <machparam.h>
  23.  
  24. #if BYTE_ORDER==BIG_ENDIAN     
  25. #define MSW 0
  26. #define LSW 1
  27. #endif
  28. #if BYTE_ORDER==LITTLE_ENDIAN
  29. #define MSW 1
  30. #define LSW 0
  31. #endif
  32.  
  33.  
  34. /*
  35.  *----------------------------------------------------------------------
  36.  *
  37.  * isnan --
  38.  *
  39.  *    Return whether a double is equivalent to NaN.
  40.  *
  41.  * Results:
  42.  *    1 if the number is NaN, 0 otherwise.
  43.  *
  44.  * Side effects:
  45.  *    None.
  46.  *
  47.  *----------------------------------------------------------------------
  48.  */
  49.  
  50. int
  51. isnan(value)
  52.     double value;
  53. {
  54.     union {
  55.     double d;
  56.     long l[2];
  57.     } u;
  58.  
  59.     /*
  60.      * Put the value into a union so we can check out the bits.
  61.      */
  62.     u.d = value;
  63.  
  64.  
  65.     /*
  66.      * An IEEE Std 754 double precision floating point number
  67.      * has the following format:
  68.      *
  69.      *      1  bit       -- sign of Mantissa
  70.      *      11 bits      -- exponent
  71.      *      52 bits      -- Mantissa
  72.      *
  73.      * If the exponent has all bits set, the value is not a 
  74.      * real number.
  75.      *
  76.      * If the Mantissa is zero then the value is infinity, which
  77.      * is the result of division by zero, or overflow.
  78.      *
  79.      * If the Mantissa is non-zero the value is not a number (NaN).
  80.      * NaN can be generated by dividing zero by itself, taking the
  81.      * logarithm of a negative number, etc.
  82.      */
  83.  
  84.     /*
  85.      * check the exponent
  86.      */
  87.     if ((u.l[MSW] & 0x7ff00000) == 0x7ff00000) {
  88.  
  89.     /*
  90.      * See if the Mantissa is zero.
  91.      */
  92.     if ((u.l[MSW] & ~0xfff00000) == 0 && u.l[LSW] == 0) {
  93.         /*
  94.          * Infinity.
  95.          */
  96.         return (0);
  97.     } else {
  98.         /*
  99.          * NaN.
  100.          */
  101.         return (1);
  102.     }
  103.     } else {
  104.     /*
  105.      * Normal.
  106.      */
  107.     return (0);
  108.     }
  109. }
  110.  
  111.  
  112.